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SUMMARY 

A description and an analysis of the thrust detection systems aboard the SERT I 
(Space Electric Rocket Test) ion thrustor flight test are given. Two types of thrust de- 
tection were employed: (1) sun sensor spacecraft spin period detectors and (2) a radial 
accelerometer centripetal acceleration detector. The thrust of a 22-millinewton ion 
thrustor produced a rate of change of spacecraft spin of approximately 0. 01 percent per 
second. This change was measured by both detection systems to yield direct thrust 
measurements agreeing to within 5 percent for 10- to 30-second averages. The radial 

accelerometer was also used as a very sensitive precession detector. The precession 
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half-cone angle of the spacecraft was measured to be 5. 3x10 radian (0. 03 ). 


INTRODUCTION 

The flight of SERT I, directed by the Lewis Research Center, on July 20, 1964, 
marked the first successful production of thrust by an ion thrustor in space. This report 
discusses the instrument systems used to measure the thrust produced. Included will be 
a description of the systems, an analysis of their function and accuracy, and a summary 
of the results obtained from the data produced by them. A report of the overall results 
of the SERT I flight has been published (ref. 1); thus, only the features of the experiment 
salient to the subject at hand will be discussed. 

The primary objective of the SERT I flight experiment (refs. 1 and 2) was to confirm 
that an ion beam could be neutralized in space. A contact- ionization thrustor using 
cesium propellant and an electron-bombardment thrustor using mercury propellant, 
which produced maximum expected thrusts of 7 and 25 millinewtons, respectively, were 
to be the first ion thrustors operated in space. In order that beam neutralization might 
be demonstrated most conclusively, the spacecraft was designed to permit direct mea- 
surement of thrust. The presence of thrust for an appreciable time would be impossible 
were neutralization lacking. 


Prior to the SERT I flight some indirect measurements of ion thrustor thrust had 
been made during vacuum-chamber testing. The method used was to place a pendulum- 
type target in the ion stream and detect the deflection produced by the impinging beam 
(ref. 3). Agreement was within 10 percent between the measured thrust values and those 
theoretically calculated on the basis of electrical measurements. 

The spacecraft was launched into a suborbital trajectory providing approximately 
1 hour of experiment time. Vacuum-chamber testing of the thrustors indicated that high- 
voltage breakdown would probably occur and cause the thrust to be intermittent, with 
bursts as short as a few seconds duration. This fact plus the low expected thrust level 
were the conditions which established the basic configuration of the SERT I spacecraft 
such that a measurement of thrust, hence beam neutralization, could be made. 

The two ion thrustors were mounted so that their thrust would cause a torque about 
the spin axis of the spin- stabilized spacecraft, and hence, a change in spin rate. This 
would allow the detection of short low-level pulses of thrust by measuring the integrated 
effect on spin rate. Thrust was thus determined by measuring spin rate and its rate of 
change. 

There were two spin- rate detection systems employed. The first consisted of two 
independent, redundant, sun sensor, spin period detectors, each of which transmitted 
one pulse per spacecraft revolution by way of separate telemetry links to the ground 
stations. The period between successive pulses from either sensor was then measured 
with ground based precision electronic counters. The second system consisted of an ac- 
celerometer mounted so that it sensed the centripetal acceleration due to spacecraft spin. 
Its output signal was transmitted redundantly through both telemetry links. 

Two ground receiving stations were used to acquire the data, Wallops Island, 

Virginia, and Bermuda. Identical thrust measurement ground equipment was located at 
each station. 


SYSTEM DESCRIPTIONS 
Simplified Spacecraft Configuration 

The basic, simplified configuration of the SERT payload was that of a right circular 
disk spinning about its longitudinal axis (fig. 1). While this configuration differs con- 
siderably from that of the actual payload (see figs. 2 and 3), it serves as a model ade- 
quate for the thrust detection analysis. Where deviations from this model are significant 
to the analysis, corrections or allowances will be made, and these will be made evident. 

The two ion thrustors were mounted opposite each other on arms extending out from 
the disk on the x-axis (fig. 1). (Note that the axes shown in fig. 1 are fixed in the disk. ) 
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Figure 4. - Sun sensor system. 


During launch the thrustors were folded 
down (fig. 2) under the disk to fit in- 
side the nose cone. The significant 
sequence of events, as planned, was 

(1) Launch, by an NASA Scout 
rocket, through third- stage 
burnout 

(2) Nose cone ejection 

(3) Spin up of fourth stage and 
payload (with small solid fuel 
rockets) 

(4) Fourth stage burn 

(5) Fourth-stage separation from 
payload 

(6) Thrustor deployment 

(7) Contact ion thrustor burn (ap- 
proximately 20 min) 

(8) Contact ion thrustor shutdown 

(9) Bombardment ion thrustor 
burn (approximately 20 min) 

(10) Reentry 

This report will be concerned pri- 
marily with the sequence after step (6). 
With the thrustors deployed, the pay- 
load would have a spin rate (-kg ) of 
about 10 radians per second, which 
was required for spin stability, and its 
z-axis would be precessing about its 
angular momentum vector with an in- 
stantaneous half-cone angle 6 . (Sym- 
bols are defined in appendix A. ) There 
were energy dissipating precession 
dampers on the payload to reduce 9 
to a maximum of 0. 035 radian (2°). 

For smaller 8 the dampers were in- 
operative. 

The thrust vectors of the two ion 
thrustors were both parallel to the y~ 
axis and in the +y-direction, and the 
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result would have been a -z-torque (i. e. , despin) from the contact ion thrustor and a +z- 
torque (i. e. , spinup) from the bombardment ion thrustor. There would be a despin fol- 
lowed by a spinup. The basic problem was to determine the thrust of the thrustors by 

measuring the very slight rate of change of co . 

z 

The following sections will describe the systems used to make this measurement. 
First will be described the sun sensor and its associated system and then the accelerom- 
eter with its associated system. 


Sun Sensor 

Sun sensor system . - There were two sun sensors located approximately ir radians 
apart on the spacecraft base plate. Each sent its signal through separate flight equip- 
ment and a separate telemetry system. At each ground receiving location either sun 
sensor signal could be selected, the choice resting upon signal quality. Figure 4 shows 
a block diagram of one of the sun sensor systems. The flight components, separate but 
identical for each sensor, consisted of a sun sensor, pulse shaper, subcarrier oscilla- 
tor, transmitter, and antenna. The ground based equipment at each receiving station 
consisted of an antenna, a receiver, a discriminator, two precision electronic time- 
interval meters, and a printer. 

The sun sensor, which produced 1 pulse per revolution of the spacecraft with respect 
to the sun, was mounted on the periphery of the spacecraft base plate, as shown in fig- 
ures 1 to 3. As the spacecraft revolved about its spin axis, sun rays passed through 
the narrow vertical slit (fig. 5) in the front of the sensor housing and illuminated a photo- 
voltaic cell mounted at the rear of the housing. The output of the photovoltaic cell was 
amplified and used to trigger a monostable multivibrator, which, in turn, was used to 
modulate a standard IRIG voltage -controlled subcarrier oscillator. This subcarrier was 
then one of the channels of the telemetry link. 

On the ground, after reception and discrimination, the signal was a replica of the 
monostable multivibrator output. This pulse train was sent through a dual -counter 
coupler in such a way that the period between successive pulses was measured alternately 
with two precision time-interval meters. The output of these meters and real time were 
recorded on a printer. 

Launch t ime e ffect on sun sensor. - To ensure illumination of the sensor by the sun 
and to prevent illumination by the Earth or moon, the time of launch was chosen to opti- 
mize the angle between the spacecraft spin axis and the line joining the sun and space- 
craft. As the spacecraft spun about its axis of maximum moment of inertia, the geo- 
metric line of sight limits of the sun sensor swept out a triangle of revolution (fig. 5). 

The upper limit of the triangle of revolution was determined by the length and location of 
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Figure 5. - Sun sensor. 



Figure 6. - Simplified accelerometer. 
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the slit aperture with respect to the photocell at the rear of the sun sensor housing. To 
minimize thrust error due to nutation and precession, the angle between the sun and spin 
axis should have been tt/2 radians, but the trajectory and, hence, the attitude of the 
spacecraft with respect to the Earth was optimized for maximum experiment time. The 
3a limit on the angle between the spin axis and a line from the photocell through the slit 
to the Earth’s horizon fixed the lower limit of the triangle of revolution at 1. 5 radians 
(85°). The upper limit of the triangle of revolution was fixed at 7r/3 radians (60°) by 
the housing geometry. The angle between the sun and the spin axis was computed as a 
function of Greenwich mean time (see appendix B). The launch time was chosen to place 
the sun as near as possible to 1. 5 radians (85°) and to place the moon outside the tri- 
angle of revolution. 


Accelerometer System 

Accelerometer. - The accelerometer was one of the type known as "vibrating 
string” instruments. A schematic view of this type of instrument, in which the strings 
are replaced by ribbons, is shown in figure 6. _ 

If an acceleration a acts along the sensitive axis A of the instrument, then a ten- 
sion force Fg, equal to ma, will be exerted on the thin metallic ribbon (m = m^ or 
m 2 ). The mechanical resonant frequency of a string in tension is 

t r mB ^a rg i 

where B is a constant and Fj is an initial tension supplied by the spring. Thus, the 
resonant frequency of the string is a function of the acceleration. If two of these string- 
restrained mass systems are opposed (see fig. 6), then the tension in one increases 
while the tension in the other decreases for a given direction of applied acceleration. 
Thus, 

f r, 1 = B 1 V F I + F S 

f r, 2 = B 2 V F I " F S 

where 

f r j resonant frequency of ribbon attached to m^ 

f r 2 resonant frequency of ribbon attached to m 2 

Bj,B 2 constants 
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As shown in figure 6, each ribbon is made flat and placed in a magnetic field in such 
an orientation that electric current passed through it will cause deflection perpendicular 
to the field. Each ribbon-magnet system is used as the tuning element in an oscillator 
circuit whose output is a sinusoidal signal with frequency a function of acceleration. The 
signals from the two oscillators are mixed and filtered to pass just the difference fre- 
quency given by 


f r, 1 " f r, 2 ” B 1 V^I + F S " B 2 V F I " F S 
Careful control during manufacture yields 

B i = b 2 = B 


So that application of the binomial expansion to equation (1) gives 

f = f i - f 2 = B (fJ/ 2 + - Ff 1//2 F S - I F- 3 / 2 F 2 + -1 Ff 5 / 2 F 3 - . . . + 

r, 1 r,2 II 2 I S g I S lg I S 


f 1/2 1 f -1/2 f 1 f - 3/2 f 2 1 f - 5/2 f 3 

-Fj +-Fj F S + g F I F S + ^ F I F S 


. + 


F I 


With Fj » Fg for the design span of the instrument 


f = B 


Pi 


( 2 ) 


or the output frequency is directly proportional to the force caused by acceleration. 

Electrical system. - One of the restraints on the electrical system, imposed by the 
payload design, was that standard IRIG FM/FM telemetry was to be used. If an accelero- 
meter had been selected whose output was a dc voltage level proportional to acceleration, 
then this signal would have been used as the input to a standard IRIG voltage -controlled 
oscillator. Since the very slight rate of change of this signal would have been swamped 
by the drift of a voltage -controlled oscillator, this method was impractical (see the sec- 
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tion Accelerometer Analysis, p. 16). This is precisely the reason for the selection of 
the particular type of accelerometer that was used. Its output, a sinusoidal voltage 
whose frequency was proportional to acceleration, was well suited to the telemetry since 
its output looked like a subcarrier channel. 

The basic quantity to be measured was the time rate of change of this frequency, 
which, as will be shown, was directly proportional to thrust. It will be shown in the 
section Accelerometer Analysis (eq. (25), p. 19) that this frequency -time function was 
modulated by a sinusoidal term whose amplitude was proportional to tan 0, where 0 is 
the precession half -cone angle. The overall or long term frequency -time function was 
then a ramp, whose slope was proportional to thrust, with a sinusoidal modulation whose 

amplitude was proportional to tan 0 (see 
fig. 7). 

The electrical system was required 
to perform three functions: 

(1) Accurately measure the slope of 
the frequency -time function ramp 

(2) Measure the amplitude and fre- 
quency of the modulation 



Figure 7. - General form of accelerometer output (not to scale). 


Flight components 



Figure 8. - Accelerometer electrical system. 
















(3) Provide a real time indication of thrust 
A block diagram of the electrical system is shown in figure 8. 

The accelerometer signal was mixed with the standard subcarrier outputs and used 
to modulate the transmitter. On the ground, the received signal was demodulated and 
passed through a constant time delay filter to separate the accelerometer signal from the 
other data channels. The sinusoidal accelerometer signal at point 1 was then reproduced 
at point 2. 

This signal was applied to a pulse shaper whose output was a short rise time pulse 
occurring at the time of each positive -going zero crossing of the input. This pulse was 
one of constant amplitude and width, so that the time average value of the pulse train 
voltage was directly proportional to frequency. 

Accurate digital measurement of ramp slope: An accurate measure of the ramp slope 
was achieved by using a preset counter to put out 1 pulse for every 400 that were put 
in so that the time -interval meters could digitally measure the total period of time (in 
psec) required for 400 cycles of the accelerometer signal to occur. Since the expected 
frequency was about 200 hertz, this yielded an average of the frequency over approxi- 
mately 2 seconds. The two time interval meters were used in alternation in order to 
avoid the loss of every other 400- cycle interval. 

These time interval data were then printed along with total time from the total time j 
meter. The time rate of change of period could then be obtained from these data and, 
when corrected for precession effects, if necessary, thrust could be calculated. 

Precession measurement: A plot of the modulation on the frequency-time signal was 
obtained by sending the output of the pulse shaper through a band-pass filter (0. 2 to 3 Hz) 
and thence to a strip chart recorder. Both modulation amplitude and frequency could be 
obtained from this. In addition, each time the time-interval meters completed a reading, 
a time synchronizing pulse was sent to another channel on this same recorder so that the 
digital data could be corrected for the effects of precession. This would be necessary, 
if there were appreciable precession, because the time-interval meter sampling times 
would not always correspond to equivalent levels on the modulation, or, in other words, 
aliasing would be introduced. 

Real time thrust indication: It was also required that an indication of thrust be avail- 
able during flight. This was accomplished by passing the pulse train through a low -pass 
filter (0 to 0. 1 Hz), the output of which was applied to a dc bias voltage in series with a 
high-sensitivity strip chart recorder. During the flight the slope of this record, and 
hence the thrust, was read with a protractor. 
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ANALYSIS 


The basic goal of the thrust detection system was to measure the thrust to an accu- 
racy of about ±5 percent with as short an averaging time as possible (10 to 60 sec) so 
that a plot of thrust against time could be obtained with reasonable thrust and time reso- 
lution. Since the thrust of the contact ion thrustor was expected to be 7. 1 millinewtons, 
whereas that of the bombardment ion thrustor was expected to be 25 millinewtons, the 
detection system, if sufficiently accurate to measure the contact ion thrustor thrust 
level, would certainly be sufficiently accurate for the bombardment ion thrustor. Thus, 
the analysis will be based on the thrust level of the contact ion thrustor. 


Sun Sensor Analysis 

Spin detection . - The equation for thrust in terms of sun sensor period measure- 
ments is obtained from the general equation for a rotating body 

Fx t = ± (Ico _) = — w + I — (3) 
dt z dt z dt 


With the worst case spacecraft temperature rise assumed, the term containing the 
rate of change of mass moment of inertia was calculated to be less than 0. 4 percent of 
the contact ion thrustor thrust and is, therefore, considered negligible. The equation 
may be rewritten, with F assumed constant over an interval At, 


Fx, 


Aw 

= I - 

At 


Since 


and 



<j0 


2tt 


z, M 


M 


(4) 


(5) 


where M is an integer, 
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2*(P m - Pi) 


( 6 ) 


Fx, = I 


where P^ =* P-^. = P because of the very low torque to inertia ratio, and At may be 
defined as MP, where M is the number of periods between measurements. 

The thrust may then be written as 


F = AP 
x t MP 3 


( 7 ) 


Equation (7) assumes the following: 

(1) No external torques other than the engine are acting on the spacecraft. 

(2) The rate of change of inertia due to thermal expansion of the spacecraft is negli- 
gible. 

(3) The engine is thrusting normal to its moment arm and normal to the spin axis. 

(See the section on accelerometer errors, p. 21, for consideration of this effect. ) 

(4) The spin rate o> z is pure spin about the z-axis (i. e. , no precession or nutation). 

(5) The thrust is constant over the measurement interval. 

Errors. - The errors in thrust measurement by the sun sensor system may be 
thought of as mechanical and electrical, and most of them are derivable from equation (7). 
The effects of precession (assumption (4)) are not derivable from equation (7) and will be 
treated separately. 

Mechanical: The error in thrust due to mass moment of inertia uncertainties is 
given by 


= i < 1 percent 

The uncertainty of the effective thrust vector location with respect to the spacecraft 
center of gravity gives rise to an error in F of 

M .^<1 percent 

These two mechanical errors would be relatively constant through the flight and should 
contribute a probable error of about 1 percent. 

Electrical: The electrical errors are the errors of the sun sensor system itself. 
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They are the errors associated with sensing, transmitting, receiving, and measuring the 
spacecraft period. 

The error in thrust due to the error in measuring period is given by 


dF\ = 2 dP 3 dP 
F/ p “ AP " P 


and since AP « P 


dF\ _ 2 dP 
,F / p AP 


( 8 ) 


( 9 ) 


Evaluation of equation (9) requires calculation of the expected values of P and AP. 
For the contact ion thrustor, AP may be computed from equation (7) to be 


FMP 3 x+ 

AP = - 

2ttI 


( 10 ) 


Then with the parameters listed in table I 

TABLE I. - PARAMETERS 


Parameter 

Design value 

Actual value 

Mass moment of inertia 


10. 0 (kg)(m 2 ) 

about x-axis, 1 ^ 



Mass moment of inertia 


11. 3 (kg)(m 2 ) 

about y-axis, I 



y 

Mass moment of inertia 

14. 4 (kg)(m 2 ) 

14. 3 (kg)(m 2 ) 

about z-axis, I 



k 

0.3 


Thrustor moment arm, x^ 

0. 557 m 

0. 557 m 

Sun sensor period, P 

0. 628 sec 

0. 7048 sec 

Initial angular velocity 

10 rad /sec 

8. 97 rad/sec 

about z-axis, u> z 



Initial half -cone angle, 8 

3. 5X10' 2 rad (2°) 

5. 3xl0" 4 rad (0. 03°) 

Thrust of contact ion 

7. 1 mN 


thrustor, F c 



Thrust of bombardment 

25 mN 

22 mN 

ion thrustor, F^ 



Accelerometer location, R 1 

0. 262 m 

0. 262 m 

Accelerometer scale factor, K 

7. 83 (H z )(sec 2 )/m 

7. 83 (H z )(sec 2 )/m 

Cycles of accelerometer signal, N 

400 

400 
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AP = 10M (jusec) 


( 11 ) 


The variation between successive period measurements dP was determined experi- 
mentally. A pulsed light source was used to energize the sun sensor at a rate of 2 pulses 
per second. The voltage pulse that fired the light was connected to the start input of the 
time interval meter system. The stop input was energized by way of light pulse, sun 
sensor, telemetry link, and discriminator output. Essentially, the time delay of the sun 
sensor system (a few msec) was being measured. Data were taken over periods of 
minutes with the spacecraft in the vacuum chamber with and without ion engines thrusting. 
Plots of the data for any given run show a maximum spread in the time delay of ±25 micro- 
seconds. However, close analysis of the data showed that the majority of successive 
readings differed by only about ±5 microseconds, and a much lower frequency drift was 
causing the wider spread in data. 

When a worst case, dP of ±5 microseconds, is assumed, the error in thrust is given 
from equation (9) to be 


/dF\ = 2X5_ .J. (i 2) 

\ F/ p 10M M 

which says that the thrust must integrate over 10 periods or 6. 3 seconds to reduce the 
error to ±0. 7 millinewton or 10 percent of the full-scale thrust. 

Precession: The assumption that the spacecraft is not precessing or nutating leads 
to two possible errors. The first would be an apparent scatter in the period measure- 
ment introduced by precession. This scatter should be periodic with time and, there- 
fore, removable from the data. The second arises from the fact that, if the spacecraft 
is coning initially, and if the cone angle damps out during the thrust measurement, then 
the energy contained in the precession mode may couple into the spin axis and cause an 
appar ent thru st . 

An analysis of this second source of error begins by assuming cylindrical symmetry 
so that the kinetic energy of rotation is 

KE ' J Ia, z + J I n(“x + w y) < 13 > 

If all of the energy in the precession mode would couple into the spin mode and there . 
were no dissipation, then, for 

0^ = initial cone angle 
Sg = final cone angle 
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I 


and 


J = J cos 6 
z 


the kinetic energy would be 


KE -j = ~ Io> 2 ifl + -LtanV] 
2 I N 


KE„ = i I w 2 0 (l +— tan 2 0„ | 

2 *. 2 l i N 


With 62 = 0 equations (15) and (16) reduce to 


2 I 2 2 

w z,l— tan 0 1 = w z,2 ' "z,l 
a N 


Simplifying and assuming w z 1 = w z 2 and tan 0 f = = 0. 035 radian give 


(14) 


(15) 


(16) 


(17) 


w z T 2 
Au)„ - -5 J- 

z p t r 


(18) 


= 1° (1. 34)(0. 035) 2 
2 

= 83x10”^ rad/sec (19) 

If the coupling took place uniformly over 500 seconds (about one -half the expected oper- 
ating time of each thrustor), the angular acceleration would be 

AC(J_ c n 

= 1. 7x10" 0 rad/ sec z 

At 

-4 2 

as compared with 2. 5x10 radian per second for the contact ion thrustor and 

-4 z 2 

8.8x10 radian per second for the bombardment ion thrustor. The error in F would 
be 
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7 percent 


( 20 ) 



1. 7xl0~ 5 

2. 5xl0' 4 


for the contact ion thrustor or 2 percent for the bombardment ion thrustor. 


TABLE D. - SUMMARY OF CALCULATED ERRORS 
IN SUN SENSOR SYSTEM 


Source 

Error 


In thrust of 

In thrust of 


contact ion 

bombardment ion 


thrustor, 

thrustor, 


F c 

F b 

Mass moment of inertia 

< 1 percent 

< 1 percent 

about z-axis, I z 



Thrustor moment arm, x. 

< 1 percent 

< 1 percent 

a ** 

Sun sensor period, P 

±.7 mN 

±. 2 mN 

Half-cone angle, 0 b 

7 percent 

2 percent 

Ultimate error with long 

1 percent 

1 percent 

averages 




a dP = ±5 jisec; M = 10. 

^Damping from 0. 035 rad to 0 in 500 sec. 


Summary: The error in thrust 
measurement for the sun sensor system 
depended chiefly on the scatter of the 
period measurement and the possible 
coupling of energy from precession to 
spin (table n). If the scatter were 
considerably less than the measured 
±5 microseconds or if longer averages 
could be taken, this source of error 
would be reduced. If the rate of pre- 
cession energy transfer were consider- 
ably less than the assumed value, this 
error would be reduced. The minimum 
error expected was that contributed by 
the mechanical parameters or ±1 per- 
cent. Even with no energy -transfer 
effects and period-measurement 
scatter, this ±1 percent error would be 
present. 


Accelerometer Analysis 

An additional thrust detection system was incorporated both as a backup for the sun 
sensor system and as a means of precession detection. The single-axis accelerometer 
described previously was selected to provide these functions. 

Location and orientation. - The accelerometer was located at point R (fig. 1, p. 3) 
with position vector R. With the accelerometer at R it was necessary to determine the 
optimum orientation of the sensitive axis with respect to R and the optimum location of 
R. The first case considered was y„ = R and x = z 0 ^ 0 with the sensitive axis per- 
pendicular to R and parallel to the x, y-plane. In this case, the accelerometer would 
have been used to sense R u> z , which is a direct measure of thrust (F = Iu> z /Xj.). How- 
ever, this orientation is impractical because of the tolerance on alinement of the sensi- 
tive axis and of the thrust axis. There would be a cross-axis centripetal acceleration on 
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the accelerometer equal to Roj z - Even if the instrument had zero cross-axis sensitivity, 
a misalinement of the sensitive axis from perpendicularity to R of only 5 microradians 
(1 arc sec) would result in a component of this centripetal acceleration along the sensitive 

— fi 9 

axis equal to 5x10 ° Ra> . Taking the ratio between this undesired component and the de- 

z 2 

sired Rcb component gives 5x10 a > /Co . Since u> ^10 radians per second and 
z 4 9 z z z 

cl> z — 2. 5x10 radians per second' 3 for the contact ion thrustor, this ratio is 2, which 
means that the undesired component would be about twice as large as the desired com- 
ponent. This is the result even under the most optimistic assumptions of zero cross- 
axis sensitivity and axis alinement to within 5 microradians. 

2 

The alternative was to measure Ra>„, the centripetal component, directly and from 
this detect the rate of change which would be proportional to thrust. The requirements 
in this case were for an extremely stable instrument since it was to detect accurately a 
very small change in a very large quantity. Taking the ratio of the change in accelera- 
tion to the acceleration for an averaging time At of 10 seconds shows this clearly: 

2Rc o_cb_ At 2d) At 

z z z 

RwJ w z 


_ 2 x2. 5xlQ- 4 x_10 _ 5 0x1Q -4 
10 

= 0. 05 percent 

Thus, the quantity of interest is only about 0. 05 percent of the quantity being measured 

directly. Since the requirement was for 5 percent accuracy on thrust, the required 

-4 5 

accuracy in the measurement of a change in acceleration was 0. 05x5x10 = 2. 5x10 = 

0. 0025 percent. The accelerometer used for this was described in the section Acceler- 
ometer (p. 7). It was sufficiently stable to meet these stringent requirements. 

Once it had been established that the measurement of cl> by means of detecting the 
rate of change of centripetal acceleration was feasible, it remained to determine the 
optimum orientation of the sensitive axis with respect to R because R might not lie in 
the x, y -plane. The sensitive axis could be oriented either parallel to R or perpendicular 
to the z-axis. The orientation chosen should be that which would enable more direct 
measurement of the centripetal acceleration and produce the more readily interpretable 
function of thrust. 

Consider the case where the sensitive axis is parallel to R. The radial acceleration 
at R would be 
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( 21 ) 


a R " 


R • R 
R 


R 


R = wXR+idX( w XR) 

_ cl) • (R x R) [o)(R • o>) - R(o) • a?)] • R 


R 


R 


( 22 ) 

(23) 


2 2 / ^ 

a R = -Ro) sin Z_(o>, R) 


(24) 


The significant feature of this result is that a R would depend only on co and not on cl>. 
This fact is important in that it would make the accelerometer signal more readily in- 
terpretable. The time rate of change of a R would be directly proportional to thrust for 
small change in u> and not complicated by dependence on co. 

If the sensitive axis were alined other than along R, an equation for the\acceleration 
along the sensitive axis, similar to equation (22), would result. However, the- x term pro- 
portional to co would not be zero because there would be no R x R multiplier as in 
equation (22). The measured acceleration would be a function of co, and the time deriva- 
tive of this acceleration would be a function of co. Thus, the interpretation of this accel- 
eration function in terms of thrust would be more complex. This analysis led to the 
orientation of the accelerometer parallel to R. As can be seen from equation (23), the 
ideal location of R would have been in the x, y-plane since this would maximize 
sin u co, R). However, physical restraints of the payload prevented this exact location; 
that is, there was insufficient room on the plane so that the accelerometer had to be 
mounted under the disk. 

Acceleration as func tion of time. - With the accelerometer located at R and 
oriented parallel to R, the acceleration a R as a function of time along that vector at 
that point had to be determined. 

The basic equations used to determine a R are Euler f s dynamical equations 

L x ~ *x^x + Cz " Jy^z^y 


L y V^y " y w z w x 


L = I <jj +(I - I )co co 
z z z ' y x' x y 


The solution to these three equations for a R as a function of time was obtained on a 
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digital computer. However, because of the complexity, a few simplifying assumptions 
will be made at this point to yield an analytically obtainable and interpretable solution. 
Features of this solution that differ significantly from the computer solution will be made 
evident in the section RESULTS. 

The simplifying assumptions are 

(1) I x = Iy (cylindrical symmetry). 

(2) ci>„ = 0. Even though co_ was the basic quantity to be measured, w„ was a very 

z z z 

slowly varying function of time because of the small torque of the ion thrustor and will be 
considered constant over a few precession periods. 

(3) L x = Ly = 0. This is really not an assumption if perfect alinement of the ion 
thrustors is assumed and there are no extraneous torques. 

Appendix C shows the derivation of ap, the result of which is 


a R - Rw z 


cos 2 /3 - 1 tan g sin 2 g cos(k t _ 6) 


N 


(25) 


Equation (25) can now be written in the form 

a R = R'<u 2 [l - Cj cos(ku> z t - 6)J 


(26) 


where 


C 1 =2 — tan 6 tan /3 
J N 


0. 046 


(27) 


k 



- 1 


(28) 


R* = R cos 2 ;3 


(29) 


which is seen to be a steady acceleration with a ripple term whose amplitude is approxi- 
mately 5 percent of the steady state (see eq. (27)). The presence of the angle 6 is seen 
to do nothing more than shift the phase of the ripple as seen by the accelerometer. 

Thrust measurement. - In the last section it was shown that 


a 


R 



(30) 
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if the small ripple term is neglected. In actuality this term could not be neglected, but 
an analysis based on its being zero was used and the result was corrected for the ripple. 

The instrument used was a single -axis accelerometer whose output was a sinusoidal 
electrical signal whose frequency was a direct, linear function of acceleration (see 
eq. (2)): 


f = Ka R (31) 

Deviations from this response were shown to be negligible in the section Accelerometer 
(P. 7). 

The equation relating this frequency to thrust can then be obtained: 


F = — (32) 


where a = d> z 

f = KR’u)5 

Z 


f = 2KR’a> z a (33) 


F = (34) 

2x^KR <x> z 


F = (35) 

2x t VKR'f 


With 


K = 7. 83 ( H z K sec2 ) 
m 


(36) 


R’ = 0.262 m (37) 

the initial frequency would be 

f Q 205 Hz (38) 
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With 


F = 7. 1 mN 

V 

(39) 

I = 14. 4(kg)(m 2 ) 

(40) 

x^ = 0. 557 m 

(41) 


the angular acceleration resulting from full-scale contact ion thrustor thrust would be 

a = 2. 75xl0~ 4 rad/sec 2 
so that for a 10- second interval (At = 10) 

j At = 2a -^ = 5. 5xl0 -4 


(42) 


(43) 


The accuracy requirement, or 5 percent of this, is about 0. 0028 percent. Thus, it has 
been shown that the accuracy required for measurement of f was 0. 0028 percent for 
10-second averaging time to 0. 015 percent for 60-second averaging time. 

Errors . - The errors in thrust measurement by the accelerometer system can be 
divided into two categories, mechanical and electrical. 

Mechanical: In order to take into account the inaccuracies of the locations of the 
accelerometer and thrustor, equation (35) can be rewritten 

F = K (44) 

2xj. sin ^_(F,x^) VlCtR' cos /_(A, R') 

where A is the sensitive axis of the accelerometer. Then 



di dXj d /_(F, xj.) 

I x t tan [_{¥, x t ) 


1 dR* | 1 

2 R’ 2 


[tan [_(a, R r )]d ^(A, R') 


It was determined by prelaunch measurements that 


di^^t^dR: <0 01 
I Xt R’ - 


(45) 


( 46 ) 
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For small misalinements the remaining two terms involving angular misalinement can be 
rewritten 


d^(F,x t ) 

tan 




(47) 


tan /_(F, x t ) 




(48) 


- [tan /_(A, R f )Jd /_(A, R') = i [/JX, R f )] : 
2 2 


(49) 


These terms would contribute only 0. 5 and 0. 25 percent, respectively, for misalinements 
of 0. 07 radian. The actual alinement was assumed to be within this limit for the thrust 
vector and was known to be within 0. 02 radian for the accelerometer. 

In summary, then, the probable error due to mechanical measurement inaccuracies 
was less than 2 percent. 

Electrical: The terms in equation (44) that were ignored in the previous section are 

/dF\ _ df _ 1 df _ 

\F/ e f 2 f 

Their size can be determined as follows. 

(1) Accelerometer scale factor K. The effect 
really twofold, the term shown and the effect on f 
latter will be examined first: 

f = KR’wJ 

f =2KRVa+R'w'K (52) 

Z z 

It will be shown that the second term in equation (52) is negligible (less than 1 per- 
cent) compared with the first; therefore, the effect of K can be neglected. The maximum 
allowable value for K can be obtained by taking the ratio of the two terms: 


1. dK 
2 K 


(50) 


on F of an uncertainty in K is 
of a rate of change in K, K. The 


(51) 
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( 53 ) 


R t A co K 

= - 5 — < 0.01 

2KR'co a 2Ka 

Za 


K < o oi 2a M 10~ 2 x2x2. 75x1 0~ 4 
K ’ co z 10 


— < 5. 5xl0“ 7 sec -1 
K 


(54) 


(55) 


This was true for the accelerometer used and furthermore was the primary reason for 
the choice of such a high quality accelerometer. Since the accelerometer was of such 
quality, the effect of an uncertainty in K could be neglected. 

(2) Time interval meter count uncertainty dT. The remaining error in electrical 
measurement was then traceable to the measurement of f: 


'dF \ = df _ 1 df 
vF J e f 2 f 


Now 


f = 



At 


(56) 


(57) 


Thus 


df = d ( f 2 ~ f l^ _ d At 
f f 2 - f l At 


(58) 


Since the measurements of f 2 and f^ 


are separate but subject to the same errors, 


df _ 2 df _ d At 
f Af At 


(59) 


Since At was of the order of 2 seconds and d At of the order of 10~° second, the 
d At/ At term could easily be neglected: 



(60) 
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However, Af « f, and thus 



( 61 ) 


Since the time interval meters did not actually measure f but T, the period of N 
cycles of f (N = 400), it is necessary to express equation (61) in terms of T and dT: 

f = ** (62) 

T 



f dT 
T 


(63) 


Af = f At 


Af = 2 KR t o) a At 
z 


Af = 


2f a At 


an 


2 df 
Af 


anN dT 

/a 

fa AtT 2 


Since the data were being taken continuously, 

At = T 


2 df 
Af 


anN dT 
z 

faT 3 


2 df 
Af 


u>z f 2 


— dT 


“ N 2 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


Equation (70) brings to light the need for the preset counter since it can be seen that a 
fairly large N is required to bring this error to the level of about 1 percent. Neverthe- 
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less N must be kept as small as possible to allow time resolution of thrust. 

When values are applied to equation (70), 

2 df m _ 10 200 2 x 1q -6 

M 2.75X10" 4 400 2 

2 df 

0. 01 or 1 percent 
Af 

where dT was taken as the 1 count (1 /isec) uncertainty of the time interval meters. 

(3) Electrical noise. In this section the required signal -to -noise ratio of the signal 
input to the pulse shaper will be established on the basis of a maximum dT of 1 micro- 
second. In the previous section dT was taken as the 1 count uncertainty of the digital 
counters. This would be the error limit if the signal input to the pulse shaper were 
perfectly noise free. 

If the signal is assumed to be a pure sine wave (of frequency f) with additive noise, 
then it can be readily shown that the following relation exists among the voltage signal- 
to-noise ratio S/J/, the error in time of the pulse shaper trigger r, and the fre- 
quency f : 


(71) 

(72) 


S _ 1 
irfr 


(73) 


With f c* 200 

- fi 

r = 10 (the resolution of the counters) 

— = 16xl0 2 or 64 db (74) 

Equation (74) shows that the required signal -to -noise ratio was about 64 decibels. 

(4) Variation in filter time delay dr^. It will be shown in this section that the time 
delay of the low-pass filter (see fig. 8, p.9) had to be constant (within a tolerance) over 
the frequency range of interest. 

If the time of the first positive zero crossing of the signal at the input to the filter is 
taken as zero, the result is a positive zero crossing at the output at 


t 


i ■ T f, i 


(75) 
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where ^ is the time delay of the filter at frequency f j. 

9 fh 

The N positive -going zero crossing emerges from the filter N cycles later at 

time 


t 2 =T + T f,2 


3tj 

T f o = Tf 1 + Af 

x > 1 3f 


3r f 

t 9 - = T + — Af 

2 1 3f 


(76) 

(77) 

(78) 


From equations (65), (68), and (78) 


*2 - h = 


2KR’c o„a 


(79) 


In order that tg - tj be considered equal to T within 1 microsecond, 

TABLE IE. - SUMMARY OF CALCULATED ERRORS 
IN ACCELEROMETER SYSTEM 


Source 

Error 

Mass moment of inertia 

In thrust of 
contact ion 
thrustor, 

F c> 

percent 

<1 

In thrust of 
bombardment ion 
thrustor, 

F b> 

percent 

<1 

about z-axis, I 



Thrustor moment arm, x^. 

<1 

<1 

Accelerometer location, R' 

<1 

<1 

Thrustor aline ment, / ( F. x t ) 

<.5 

<. 5 

Accelerometer aline ment Z_ (A, R f ) 

<. 25 

<.25 

Accelerometer scale factor, 

<1 

<1 

K, and K 



Frequency of accelerometer signal, 

<1 

<1 

f, or period of N accelerometer 
signal cycles, T a 



Ultimate error expected 

<5 

1 


a For dT = ±1 jxsec; N = 400; S/jf = 64 db. 
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( 80 ) 


dTj; 

— -< 45 /J.sec/Hz 
df 

The filter was a special design which met this requirement. 

Summary: The total probable error in thrust measurement for the accelerometer 
system was expected to be 3 to 5 percent of the thrust of the contact ion thrustor 
(table m). 


RESULTS 

The SERT I spacecraft was launched July 20, 1964, in a ballistic trajectory. The 
contact ion thrustor failed to start because of a high-voltage breakdown. The bombard- 
ment ion thrustor operated for a total of about 30 minutes, during which time it was 
stopped and successfully restarted a number of times (ref. 1). While operating it pro- 
duced a measured thrust of about 22 millinewtons . 

The fact that only the larger of the thrustor s operated made the task of thrust detec- 
tion less critical than would have been the case if the smaller thrustor had functioned. 

The following sections will discuss the results obtained from the thrust detection systems. 


Sun Sensor Data 

Figure 9 is a plot of spin period against time as determined by the sun sensor 
shortly after arm deployment. These data are typical of those obtained during periods of 
zero thrust throughout the flight. During periods of thrust the average value of period 
sloped markedly upward. The sun sensor period oscillated about its average value of 
0. 70480 second at 0. 47 hertz or one-third the spin frequency with a peak-to-peak devia- 
tion of approximately 300 microseconds. This modulation at one-third the spin frequency 
was caused by precession and remained throughout the flight, while the peak-to-peak 

deviation decreased almost linearly 
to about 100 microseconds at the end 
of the flight. 

Because the period data were 
modulated at one -third the spin fre- 
quency, the thrust was computed by 
letting M equal 3 in the thrust equa- 
tion (eq. (7)). The difference be- 
tween successive third periods 


.7050 


.7048 


. 7046] 



05 5:07 5:09 5:11 5:13 

Time after launch, min.-sec 

Figure 9. - Sun sensor period as function of time. 
Wallops Island ground station. 


5:15 


5.17 


Data taken from 
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showed a jitter of about 5 microseconds, in agreement with that observed during pre- 
flight ground testing. This gives an error due to jitter of ±1.7 millinewtons or ±8 percent 
of full-scale bombardment ion thrustor thrust over a measurement time of only 2 seconds. 
Because the contact ion thrustor failed to work, the spacecraft coasted for several min- 
utes. The spin rate during coast was observed to be constant within the observable 
limits, which verified the assumption of negligible external torque and indicated that 
there was negligible energy coupling into the spin axis. 

The sun sensor data provided an average thrust measurement over 2 seconds to 
within about ±9 percent of full-scale thrust. When the thrust was constant for 30 seconds 
or more, the error was determined primarily by the mechanical tolerances on the space- 
craft. Thus, these longer averages approached an accuracy of ±1 percent. 


Accelerometer Data 

Precession. - In the analysis it was assumed that the payload was cylindrically 
symmetrical. In practice this was not the case since the actual mass moments of 
inertia were 

l x = 10. 0 (kg)(m 2 ) 

Iy = 11. 3 (kg)(m 2 ) 

I z = 14. 3 (kg)(m 2 ) 



Figure 10. - Real time thrust record. 


The principal effect of this lack of 
symmetry on the accelerometer analy- 
sis was to make invalid the result that 
6 was a constant in time. In other 
words, nutation was introduced. Be- 
cause of the complexity that nutation 
introduces into the analysis, a digital 
computer solution of Euler's equations 
was used in conjunction with the data 
from the strip chart record of preces- 
sion to determine the actual value of 6 
and its variation with time. (See the 
section on the accelerometer electrical 
system, p. 8. ) During the first portion 
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of flight (w z = 8. 97 rad/sec) 6 oscillated between 3. 5x10”^ (0. 02°) and 5. 3x10""^ radian 
(0. 03 u ) and gradually decayed to about one-third of this value in agreement with sun 
sensor data. The value of d calculated from equation (25) for the first portion of the 
flight was 7. 4x1 0"^ radian (0. 042°), which indicates that the assumption of symmetry 
was warranted. The gradual decay of 6 indicated that energy from the precession mode 
may have been coupling to the spin mode at a fairly constant rate. However, the small 
value of d indicated negligible initial energy in the precession mode compared with the 
energy change of the spin mode due to thrust. Thus, correction for the precession 
energy transfer effect on the thrust data was not necessary. 

The period of the modulation of accelerometer frequency output, as determined from 
the same strip chart record and in agreement with the computer solution, was about 
2. 1 seconds as compared with 2. 05 seconds calculated from equation (25). 

Thrust (re al time readout). - Thrust readout was accomplished in real time as 
described in the section on the accelerometer electrical system (p. 8). A section of 
the strip chart record is reproduced in figure 10. The horizontal axis is time, and the 
vertical axis is spin rate, though it should be noted that the vertical axis on this plot 
represents only a small portion of the total spin rate since most of the signal was elec- 
trically biased out. The slope of the plot is proportional to thrust. It can be seen from 
this plot that the 22-millinewton thrust was detected with ample resolution for real time 
readout. The small oscillation on this record represents that portion of the precession- 
caused modulation that the low -pass filter failed to exclude. Its period is 2. 1 seconds, 
as explained in the previous section. 

Thrust (digital readout). - Thrust was calculated as a function of time from the 
digital data recorded from the time-interval meters. Because of the very small cone 
angles that were measured, the digital data were used without correction for precession 
energy transfer. Each value of the time interval required for 400 cycles of accelero- 
meter output frequency was subtracted from the previous value. From this difference 
was calculated f for the mean point in time. Equation (35) was then used to calculate 
F. 

The major problem encountered with these thrust data was the fact that, during 
certain portions of the flight, the sampling period (i. e. , the period of 400 cycles) was 
just a bit longer than the period of precession. This led to large aliasing error. Fig- 
ure 11 (p. 30) shows a plot of thrust against time for the period T + 27 minutes to 
T + 28 minutes, a period chosen to illustrate the large effects of aliasing. From fig- 
ure 11 it can be seen that the uncorrected accelerometer data, averaged over a minute, 
agree with the sun sensor data. A point -by -point correction for this aliasing, obtained 
by applying equation (25) to the data from the strip chart record of precession, yielded 
corrected accelerometer data (also shown in fig. 11) which agreed very well with the 
sun -sensor data. 
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Time after launch, minzsec 


Figure 11. - Accelerometer and sun sensor thrust data with large alias- 
ing effect 



Time after launch, min:sec 


Figure 12. - Accelerometer and sun sensor thrust data with negligible 
aliasing effect 
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Figure 12 shows a plot of thrust against time for the period T + 45 minutes to 
T + 46 minutes, a period chosen to illustrate the results obtained when negligible aliasing 
was observed. It can be seen that the accelerometer data and sun-sensor data agree 
during this period even through a short period of engine shutdown. The accelerometer 
data were not corrected for aliasing in figure 12. 


CONCLUDING REMARKS 

From the data presented in the preceding sections it is seen that the thrust of the 
bombardment ion thrustor was measured with agreement to within 5 percent between the 
sun-sensor and accelerometer thrust detection systems. That the two systems comple- 
mented each other is evident from the fact that the accelerometer data showed the pre- 
cession angle to be sufficiently small so that the sun-sensor data could be treated as 
being relatively free from the energy transfer effect. The sun-sensor data were used, 
on the other hand, to make an in-flight check of the accelerometer operation during the 
period of no thrust. 

It should be pointed out, in conclusion, that the use of a radial accelerometer, as 
on this payload, provided a very sensitive means of detecting small precession angles. 
This was shown by both the simplified analysis and the results. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 5, 1966. 
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APPENDIX A 
SYMBOLS 


A 

sensitive axis of accelerom- 

I 


eter 


IS A,R) 

angle between vectors A 

% 


and R 


a 

acceleration 


a 

R 

acceleration along R 

*x’ l y’ l z 

B l> B 2 

constants 


F 

thrust 


/ dj\ 

\ F /I 

error in F due to errors in 
I (typical notation used 

J 

J z 


throughout) 


F b 

thrust of bombardment ion 

K 


thrustor 

KE 

F c 

thrust of contant ion thrustor 

k 

F I 

initial tension in accelerom- 


eter spring 

L x’ L y’ L : 

F s 

force on "proof mass" of 

M 


accelerometer due to ac- 


celeration 

m 

f 

frequency of accelerometer 

N 


signal 

f r 

resonant frequency of ac- 

P 


celerometer string 


f 0 

initial frequency of ac- 



celerometer signal 

R 


maximum principal mass 
moment of inertia 

principal mass moment of 
inertia normal to I 
(assuming cylindrical 
symmetry) 

mass moments of inertia 
about x, y, and z, 
respectively 

unit vectors along x, y, 
and z, respectively 

angular momentum 

angular momentum about 
z-axis 

accelerometer scale factor 

kinetic energy 



torque about x, y, and 
z, respectively 

number of periods of sun 
sensor data 

mass 

number of cycles of ac- 
celerometer signal 

sun sensor period or time 
of one spacecraft revo- 
lution 

location of accelerometer 
or magnitude of R 
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R' 

R cos p 

T 

S/jV 

signal-to-noise ratio 


T 

period of N accelerometer 
signal cycles 

T f 

CO 

t 

time 

CO j CO CO 

x, y,z 

principal axes of spacecraft 
(see fig. 1) 

x y 7 z 

v y a » z a 

location of accelerometer 

CO n 


(see fig. 1) 

x, 0 


location of either ion thrustor 

Subscripts: 

a 

angular acceleration 

e 

“c 

angular acceleration due to 

M 


contact ion thrustor thrust 

m 

/3 

defined by fig. 1 

Superscripts: 

6 

defined by fig. 1 

9 

half- cone angle (see fig. 1) 

C) 

a 

standard deviation 

(•),(••) 


time-interval meter reso- 
lution 

filter time delay at fre- 
quency f 

angular velocity 

angular velocities about 
x, y, and z, respec- 
tively 

u>x at t = 0 

electrical 

Mth period 

mechanical 

vector 

first and second time deri- 
vatives 
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APPENDIX B 


METHOD OF COMPUTING ANGLE BETWEEN A CELESTIAL BODY AND 
SPIN AXIS OF A SPIN STABILIZED SPACECRAFT 

The SERT I spacecraft has two sun sensor systems mounted on the periphery of the 
spacecraft baseplate. As the spacecraft spins about its axis of maximum moment of 
inertia, the sun sensor views or "sweeps out" a triangle of revolution, as shown in 
figure 13. The limits of the triangle of revolution are determined by the geometry of the 
sun sensor housing and may be expressed as angles measured from the spin axis. To 
ensure illumination of the sensor by the sun and to prevent illumination by the Earth and 
the moon, the angle between the spacecraft spin axis and the sun must be controlled. 

The trajectory of the SERT I spacecraft and, hence, the attitude of the spin axis was 
optimized for maximum experiment time. The launch window is then the time interval 
during which the angle between the spin axis and the sun is within the limits of the tri- 
angle of revolution swept out by the sun sensor. 

A general solution of the angle between the spin axis of a spin-stabilized space- 
craft and the sun or any other celestial body can be calculated as a function of Greenwich 
mean time. 

The following basic assumptions are made: 

(1) The spacecraft coordinates are fixed and known with respect to an inertial refer- 
ence system at rocket burnout and are constant because of spin stabilization. 

(2) The angle between the spin axis and the sun is constant throughout the flight, be- 
cause of the relatively small distance traveled in the ballistic trajectory compared to the 
distance from the sun. 

The solution is obtained by projecting the position of the sun, the spacecraft center 

of gravity, and the spacecraft spin axis 
on to the celestial sphere (as shown in 
figure 14) and solving for the angle be- 
tween the spin axis and the sun with 
spherical trigonometry. 

The sun may be located on the celes- 
tial sphere by finding the Greenwich hour 
angle or longitude and the declination or 
latitude as a function of Greenwich mean 
time in The American Ephemeris and 
Nautical Almanac (ref. 4). The space- 
craft center of gravity may be located by 
knowing its longitude and latitude at 
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booster engine burnout. Finally, the spin axis is located by knowing the pitch and 
azimuth at booster engine burnout. From figure 14, the following angles are defined: 


PVS calculated 

PVZ azimuth of spin axis 
SOP colatitude of sun 

SVZ calculated 

VOP colatitude of spacecraft center of gravity 

VOS angle between local vertical and sun 

VOZ copitch of spin axis 

VPS local hour angle 

ZOS angle between spin axis and sun 


From the law of cosines 

cos VOS = cos SOP cos VOP + sin SOP sin VOP cos VPS (Bl) 

From the law of sines 


sin SOP sin VPS 

sin PVS = = = (B2) 

sin VOS 
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From the law of cosines 


cos ZOS = cos VOS cos VOZ + sin VOS sin VOZ cos SVZ (B3) 


A digital computer was used to solve the preceding equations in 1-hour increments. 
Figure 15 shows the angle between the spin axis of the SERT I spacecraft and the sun 
(angle ZOS) plotted as a function of Greenwich mean time for July 20, 1964. The shaded 
portions indicate the maximum and minimum limits of the triangle of revolution swept 
out by the sun sensor. The intersection of these limit lines and angle ZOS define the 
time limits of the flight. It is apparent that there are two launch windows, one in the 
morning and one in the afternoon. When the angle between the spin axis and the moon is 
plotted, however, the moon is seen to cross through the afternoon window. 


1501 


ST 



/-Moon 






Launch 

window- 


02 


04 06 



&M I I 

08 10 12 14 

Greenwich mean time 


Figure 15, - Angle between SERT I spin axis and sun and spin axis and moon as function of Greenwich 
mean time for July 20, 1964. 



APPENDIX C 


DERIVATION OF MEASURED ACCELERATION AS FUNCTION OF TIME 

The derivation of a R as a function of time is accomplished by starting with Euler's 
dynamical equations 


L x *x^x + < J z ” ^y^z^y 

L y - yv - 

L z “ l z?z * «y - yv, 
and the assumptions explained in the text 


If these assumptions are applied and 


*x = I y 


<h z = ° 


L = L = 0 
x y 


l z =I 


*x ~ *y " % 


I - I 


k = 


N 


a N 


(Cl) 

(C2) 

(C3) 


the result is 


cb x + ka) y w z = 0 (C4) 

d> y - ko> x a) z = 0 (C5) 

Equations (C4) and (C5) are two simultaneous linear differential equations with constant 
coefficients. With the choice of u> y = 0 at t = 0 the solutions are 
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w x = w x, 0 cos k< V 


“y = w x, 0 sin ^ 


(C6) 

(C7) 


It is necessary to express a 0 in terms of a) and co . From equations (21) and (22) 

-tv. y x 


_[wx(ti)XR)]- R 
a R = R 


a R “ 


(R X o>) • (a) X R) 
R 


- _ (R X co)' 

a R r~ 


<C8) 


Neglecting the negative sign (which merely indicates that a^ is centripetal) and expand- 
ing yield 


*R = ^[ny a " z - z a"y> + !< z a“x ’ Yj + *< x a"y ' wf 


1 

a R = R 


x a < w y + “5 + y a < w x + w z> + z a ^x + ' ^a^y "z + Va^z + WVtyl 


(C9) 


From figure 1 



x„ = -R cos j3 cos 6 

<x 

(CIO) 


y_ = -R cos /3 sin 6 

a 

(Cll) 


z = -R sin 

a. 

(C12) 

Substituting equations (CIO), (Cll), and (C12) into equation (C9) gives 


a R - R 

2„ 2^/2 2\ 2 n . 2., 2 2y 

cos /3 cos 6(a) + a> ) + cos /3 sin 6(a)„ + co ) 
y z x z 



2 2 2 

+ sin 0(u) x + a>y) - 2(cos /3 sin |3 sin 6a)yCi> z 



2 

+ cos j3 sin /3 cos 6 a) a) + cos /3 sin 6 cos 6a> a)) 
x z x y 

(C13) 
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These expressions, for co^ and co x and for in terms of u>y, c o x , and u> z , must 
be combined with an expression for 8 in terms of co x q , which can be obtained by con- 
sidering the invariance of the angular momentum vector: 


J z = J cos 8 


Ico z =yi I 2 u,2 + i2 (u) 2 + a; 2 )cose 


(C14) 

(C15) 


tan 8 = 


■nV w 


2 2 
+ CO 

X T w y 


Ico. 


(C16) 


From equations (C6), (C7), and (C16) or from figure 1, 


tan 8 = 


^x, 0 

Ico_ 


(C17) 


Thus 


o> n = — co tan 8 
x. 0 t z 

a n 


(C18) 


and 


a R = Ra, z < 


cos 2 (3 cos 2 6 


, ,1 tan 0\ 2 . 2, . 

1 + (— ) sin kco z t 


N 


+ cos 2 0 sin 2 6 


I tan 8\ 2, f 

1 + / ] cos kco t 

T / Z 

X N 


+ sin 2 /3 (— tan 8 

Vn 


2 * * atl - (cos sin ft sin 5 sin kco z t + cos j3 sin /3 cos 5 cos kco^t 
Im V 


A N 

+ - tan - cos 2 /3 sin 5 cos 6 sin kco^t cos kco^t 

% 


t 


(C19) 


After trigonometric simplification 
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R 


<2 


r 

icos 2 /3 + — 

2 

(l tan 6\ r + g j n 2^ _ cog 2^ cog 2(kcu t - 6)1 

2 

V. 

V % / L 


sin 2/3 cos(ko) t - 6)) 

I z 

X N 


(C20) 


The first term gives the expected acceleration if 0 were zero; the second is a 
quantity oscillating at frequency ko) z /jr and amplitude 1/2 [(I/t^) tan 6 cos /3j 2 
the third is a quantity oscillating at one-half the frequency of the second and amplitude 

Rw 2 (^n) tan 6 sin 2/3 ' 

In order to interpret this equation, it is necessary to look at some of the relative 
magnitudes that result when the following numerical values are assumed: 

/3 = 0. 473 radian - actual value 

8 = 0. 84 radian - actual value 

0 = 0. 035 radian - maximum allowed by dampers 

R = 0. 33 m - actual value 


< jl> z = 10 rad/sec - assumed initial value 


I h h 

— = 1. 3 - approximate average of — and — 
% 1 X l y 


The ratio of the amplitude of the third term to that of the second (in eq. (C20)) yields 

2 sin 2/3 4 tan /3 


I 2 I 

— tan 6 cos /3 — tan 0 

hf % 
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(C21) 


so that the second term amounts to nothing more than a small amount of second harmonic 
distortion on the third term. 

Equation (C20) reduces to the following if the second term is neglected: 


a B* R “z 


2 a I tan 0 sin 2/3 . . 

cos /3 cos ( ko) t - 6] 

I V Z 

N 


(ku> z t - 6) 


(C22) 


40 



REFERENCES 


1. Cybulski, Ronald J. ; Shellhammer, Daniel M. ; Lovell, Robert R. ; Domino, 

Edward J . ; and Kotnik, Joseph T. : Results From SERT I Ion Rocket Flight Test. 
NASA TN D-2718, 1965. 

2. Gold, Harold; Rulis, Raymond J. ; Maruna, Frank A. , Jr. ; and Hawersaat, 

William H. : Description and Operation of Spacecraft in SERT I Ion Thrustor Flight 
Test, NASA TM X- 5 20 50, 1964. 

3. Kerslake, William R. ; and Pawlik, Eugene V. : Additional Studies of Screen and Ac- 

celerator Grids for Electron-Bombardment Ion Thrustors. NASA TN D-1411, 1963. 

4. Anon: The American Ephemeris and Nautical Almanac. U. S. Government Printing 

Office, 1964. 


NASA-Lajigley, 1966 


E-3125 


41 


"The aeronautical at2d space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof .” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



